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In  this  report  we  show  that  a non-linear,  adaptive  processor  called 
PHILTRE  is  useful  In  detecting  long-period  Rayleigh  waves  In  various  low 
slgnal-to-nolse  ratio  situations:  signals  burled  In  noise  at  various 

levels,  two  signals  mixed  at  various  azimuths  and  relative  amplitudes,  and 
a suite  of  visually  undetected  signals  at  LASA,  ALFA,  and  NORSAR  from  the 
Kurlls-Kamchatka  region. 

PHILTRE  lowered  the  detection  threshold  for  Rayleigh  waves  burled  In 
noise  by  about  6 dB.  It  was  able  to  separate  two  signals  If  their  azimuthal 
separation  was  greatest  than  60  degrees  and  If  at  the  same  time  the  amplitude 
of  the  second  signal  was  at  least  20  percent  of  the  amplitude  of  the  first 
signal.  It  levered  Raylelgh-wave  50Z  detection  thresholds  by  roughly  0.2  m,^ 
unit  at  the  three  long-period  arrays. 
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INTRODUCTION 

This  report  examines  a non-linear  adaptive  processor,  termed  PHILTRE 
and  described  below,  for  long-period  Rayleigh  waves  which  are  recorded  in 
various  low  S/N  ratio  situations.  One  situation  Involves  burying  signals 
in  nolSe  at  various  levels,  another  Involves  mixing  two  signals  at  various 
azimuths  and  relative  amplitudes,  and  a third  Involves  a suite  of  visually 
undetected  signals  from  the  Kurils-Kamchatka  region. 

The  time-varying,  non-linear  processor  (termed  PHILTRE)  used  to  enhance 
long-period  surface  waves  in  this  report  was  first  presented  by  Simons  (1968). 
This  technique  operates  on  all  three  traces  after  rotating  the  horizontal 
traces  to  be  radial  and  transverse  to  the  great-circle  path  from  epicenter 
to  station.  The  traces  are  first  transformed  to  the  frequency  domain.  Then 
the  frequency  components  are  weighted  according  to:  1)  apparent  azimuth  of 

approach  relative  to  the  expe.':ted  azimuth  and  2)  closeness  of  observed  phase 
lag  between  vertical  and  radial  traces  to  the  expected  90“  for  Rayleigh  waves. 
Finally  the  weighted  harmonic  components  are  inverse  transformed  to  produce 
the  output.  Evaluations  of  this  type  of  processor  have  been  made  by  Choy 
and  McCamy  (1973)  and  Nelson  and  von  Seggern  (1974)  using  digital  recordings 
from  the  VLPE  network.  Results  shown  In  these  studies  as  well  as  Simons' 
paper  were  encouraging;  however,  since  most  of  the  signals  In  these  studies 
were  visually  detectable  anyway  before  application  of  the  processor,  a 
stringent  test  of  the  processor's  capability  at  very  low  S/N  raf*os  has  not 
been  provided.  Such  a test,  presented  in  this  report.  Is  necessary  because 
the  non-linear  nature  of  the  processor  should  cause  its  efficacy  to  deteriorate 
rapidly  below  S/N  ratios  of  one. 


Simons,  R. , 1968,  PHILTRE  - A surface  wave  particle  motion  discrimination 
process.  Bull.  Seism.  Soc.  Am.,  v.  58,  p.  o29-637. 

Choy,  G.  and  K,  McCamy,  1973,  Enhancement  of  long-period  signals  by  time  varying 
adaptive  filters,  J.  Geophys.  Res.,  v.  78,  p 3505-3511. 

Nelson,  D.  and  von  Seggern,  D.  H.,  1974,  Signal  enhancement  of  LPE  data,  SDAC- 
TR-73-4,  Teledyne  Geotech,  Alexandria,  Virginia. 
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Procedure 


For  this  part  of  the  study  we  selected  one  event  and  three  noise  s^ples 
from  two  LPE  sites;  CTA  (Chi.rtars  Towers,  Australia)  and  KIP  (Klpap«,  Hawaii), 
The  coordinates  of  the  stations  and  events  and  the  times  of  the  event  s^imples 
are  shown  In  Tables  I and  II.  The  noise  samples  are  from  the  first  9 days  of 
June  1972.  The  frequency  renponae  of  the  LPE  system  Is  centered  at  approxi- 
mately 40  seconds  period.  Seismic  data  are  digitally  recorded  at  the  sites 
with  a sampling  Interval  of  one  second. 


The  Input  traces  aie  Initially  bandpass  filtered  for  periods  between 
approximately  12  and  80  seconds.  A sly.— Lilnute  section  of  the  LR  phase  having 
maximum  amplitude  Is  mixed  at  various  nagnlflcatlons  with  each  of  the  three 
nolso.  samples  at  that  sLatloc . Th«  n the  m jced  traces  are  procesi^ed  through 
the  program  FHILTRE. 


Results 


Figures  1 through  6 are  30— minute  plots  of  some  of  the  mixed  traces. 

Each  figure  shows  6 traces:  traces  1-3  are  Input  traces  to  PHILTRE,  bandpass 

rotated  into  vertical,  radial,  and  transverse  components,  and 
traces  4—6  are  PHILTRE  outputs.  Although  gains  are  equalized  between  the 
three  conq/onents  for  computation,  they  are  not  equalized  on  the  plots  wh&ro 
the  traces  ire  plotted  so  that  the  maximum  amplitudes  are  the  same  height  on 
each  trace.  The  amplifications  to  the  right  of  each  processed  trace  are 
correct  only  with  respect  to  one  another;  they  give  a measure  of  the  relative 
amplitude  of  the  horizontal  traces  to  the  vertical  traces.  If  It  were  desired 


to  calculate  M^  at  e.g.  20  secouds  from  these  plots  It  would  be  necessary  to 


multiply  the  Z component  time  series  by  a constant  equal  to  the  Inverse  of 
the  20-second  PHILTRE  coefficient.  Where  there  Is  a small  amplitude  output, 
the  output  traces  show  spikes  and  sudden  jumps. 

To  generate  Figure  7 which  shows  S/N  output  as  a function  of  S/N  Input  we 
first  Imbedded  signals  In  nol&e  at  high  S/N  values.  Hie  S/N  values  were  then 
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TABLE  II 

LPE  Station  Parame*'tjrs 


Station 

CTA 

KIP 


Coorrf-»:.ates 

20. IS  146. 3E 
21. 4N  158. OW 


1 
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defined  as  the  maximum  peak- to- trough  signal  amplitude  divided  by  six  times 
the  rms  noise  amplitude  for  the  noise  alone  in  the  first  256-second  window 
on  the  plot.  Successive  traces  were  generated  by  adding  the  signal  to  the 
noise  with  increasingly  smaller  signal  amplitudes.  The  S/N  values  for  these 
traces  were  calculated  by  scaling  the  S/N  values  in  proportion  to  the  signal 
amplitude  liince  it  would  be  impossible  to  empirically  measure  S/N  input  for 
those  cased  where  the  mixed  signal  falls  below  the  visual  threshold  on  input. 

Examination  of  Figure  7 shows  that  the  mean  processor  gain  is  lA  dB  at 

the  visual  threshold  which  is  at  a S/N  input  of  1.0,  which  corresponds 

roughly  to  an  empirical  50/.  M threshold  (von  Seggern,  1975).  The  curves 

8 

in  Figure  7 are  plotted  for  successively  lower  S/N  input  until,  in  the 
analyst's  judgment,  the  signal  cannot  be  detected  on  the  processed  traces. 

(The  filled  symbols  denote  detections  and  the  empty  symbols  denote  those  cases 
where  the  signal  was  not  detected.)  The  asterisks  plotted  in  Figure  7 are 
midway  between  the  last  detection  and  first  non-detection  and  are  the  estimated 
points  at  which  the  signal  would  be  detected  with  50Z  probability. 

The  star  in  Figure  7 gives  the  mean  of  the  50%  detection  points.  Although 
the  standard  deviation  of  the  mean  seems  quite  large,  the  result  is  in  agree- 
ment with  other  results  to  be  presented  below.  We  see  that  the  detection 
threshold  for  S/N  inputs  has  been  reduced  approximately  6 dB,  However,  at 
this  detection  threshold,  the  S/N  output  is  not  0 dB;  it  is  approximately 
5 dB,  The  5 dB  of  improvement  is  due  to  the  non-linearity  of  the  processor 
and  to  noise  pulses  in  the  signal  window.  As  we  have  seen,  if  S/N  inputs  are 
reduced  further  the  analyst  cannot  see  the  signal. 

Thus  the  processor  promises  a decrease  in  the  M detection  threshold  of 

s 

approximately  0,3  units.  This  result,  if  achievable  in  routine  applications, 
is  comparable  to  those  for  other  processors,  such  as  matched  filtering  and 
multichannel  filtering.  There  is  the  possibility  that  PHILTRE  could  be 
applied  to  the  output  of  such  processors  resulting  in  larger  additive  gain. 


von  Seggem,  D.  H.,  1975,  Final  report  on  the  analysis  of  recordings  from  the 
Very  Long  Period  Experimental  stations.  Teledyne  Geotech,  Alexandria, 
Virginia. 
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SIGNAL  SEPARATION  OF  EVENTS  RECORDED  AT  LPE  SITES 


Procedure 


In  this  section  we  examine  the  ability  of  PHILTRE  to  recover  a signal 
when  another  signal  Is  present.  As  explained  in  the  introduction,  PHILTRE 
weights  the  frequency  components  of  the  trace  according  to  deviation  from 
expected  azimuth.  From  0 to  90  degrees  from  the  expected  azimuth,  the 
frequency  components  are  weighted  according  to  the  fourth  power  of  the  cosine 
of  the  expected  azimuth  minus  the  calculated  direction  of  motion.  For  calcu- 
lated directions  between  90  and  180  degrees  from  the  expected  azimuth,  the 
frequency  components  are  multiplied  by  zero. 


To  study  philtre's  capability  In  the  presence  of  interference,  two 
signals  were  chosen  from  LPE  recordings  at  the  site  CTA  and  two  from  KIP. 
Table  III  lists  the  coordinates  of  the  events  with  which  the  signals  were 
associated  and  the  times  of  all  samples. 


A four-to-five  minute,  maximum-amplitude  section  of  the  LR  phase  for 
the  second  event  recorded  at  each  site  is  mixed  at  various  lags,  magnifica- 
tions, and  azimuthal  separations  with  the  first  event.  Then  the  mixed  traces 
are  rotated  to  the  back  azimuth  of  the  second  signal  and  processed  through 
PHILTRE  in  an  attempt  to  detect  the  second  signal. 


Results 


Figures  8 through  11  are  30-  minute  plots  of  the  four  signals  before 
mixing.  The  signals  have  been  processed  through  PHILTRE,  and  these  PHILTRE 
outputs  should  be  compared  with  examples  of  PHILTRE  outputs  of  the  mixed 
traces  shown  in  Figures  12  through  15.  Figures  16  and  17  summarize  detection 
results  for  all  the  cases  using  the  relative  amplitudes  of  the  first  and 
second  signals  versus  their  azimuthal  separation  for  lags  of  1 minute  and 
2 minutes  between  the  two  phases,  respectively.  A lag  of  2 minutes  represents 
about  one-half  overlap  of  the  phases  and  a lag  of  1 minute  is  about  three- 
quarters  overlap  of  the  LR  phases.  The  second  signal  was  defined  to  be 
detected  if  its  amplitude  was  at  least  twice  that  of  the  first  signal. 


-20- 

( 

i 


f 


f 


Figure  12.  llixing  of  two  signals  at  station  CTA,  showing  marginal  separation 
of  second  signal.  Lag=l  minute.  Azimuthal  separation  = 30°.  Amplitude  of 
second  signal  • 40%  of  first  signal. 
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AZIMUTHAL  SEPARATION 


Figure  16.  Relative  amplitude  vs.  azimuthal  separation  for  lag  = 1 minute 


At  azimuthal  separations  of  less  than  about  60  degrees,  the  second 
signal  is  not  separated  from  the  first  signal  unless  the  second  signal  was 
originally  at  least  twice  the  amplitude  of  the  flrct  signal.  The  second 
signal  is  also  not  detected  when  it  is  less  than  about  20  percent  of  the 
amplitude  of  the  first  signal  even  when  the  azimuthal  separation  is  greater 
t>^an  about  60  degrees.  Above  20  percent  of  the  amplitude  of  the  first 
signal  and  for  more  than  60  degrees  of  azimuthal  separation,  the  second 
signal  is  detected  while  the  first  signal  is  suppressed.  These  results 
demonstrate  the  usefulness  of  PHILTRF  in  those  cases  where  masking  of  signals 
occurs;  in  the  special  case  that  the  "noise"  is  actually  another  signal, 
PHILTRE  can  suppress  the  "noise"  considerably  if  its  azimuth  is  quite 

different  and  thereby  recover  signals  well  below  the  instantaneous  "noise" 
level. 


THRESHOLD  IMPROVEMENT  USING  PHILTRE 


Procedure 

Long-period  LR  beams  at  the  arrays  ALFA,  LASA,  and  NORSAR  for  events 
from  the  Kamchatka-Kurils , Seismic  Region  19,  in  February  and  March  1972  have 
been  formed  and  visually  analyzed  at  the  SDAC.  Eighty-seven  percent  (87%) 
of  the  possible  events  were  undetected  visually  on  the  beams  at  one  or  more 
arrays.  (Masking  contributed  to  this  high  rate  of  non-detection,  as  well 
as  the  signals  being  below  ambient  noise  levels.)  In  order  to  estimate  how 
much  PHILTRE  could  lower  the  detection  threshold,  all  beams  were  processed 
where  the  signal  was  masked  or  was  simply  undetected.  To  demonstrate  the 
effectiveness  of  PHILTRE  alone,  no  bandpass  filter  was  used  in  the  PHILTRE 
runs  for  these  cases.  Outputs  from  PHILTRE  were  submitted  to  the  same  analyst 
who  originally  read  the  unprocessed  beams,  and  he  decided  whether  the  LR 
phase  was  not  detected.  In  addition,  all  of  the  NORSAR  data  was  simply  band- 
pass filtered,  using  frequency  cutoffs  of  .03  and  .06  hertz,  and  showed  no 
new  detections;  thus  we  concluded  that  any  improvement  by  PHILTRE  over  raw 
traces  would  also  be  an  improvement  relative  to  band-pass  filtering  of  the 
raw  traces. 

Results 

Figures  18,  19,  and  20  are  samples  of  filtered  outputs.  Figures  21,  22, 
and  23  show  the  percentage  of  detection  of  LR  versus  bodywave  magnitude  for 
ALFA,  LASA,  and  NORSAR  data  individually.  Crosses  denote  the  detection  rates 
from  visual  beam  analysis,  and  circles  represent  the  rates  after  the  beams 
having  no  detection  or  masking  were  processed  through  PHILTRE  and  rechecked 
by  the  same  analyst.  Normal  distribution  curves  were  fitted  by  computer 
f ^ through  each  set  of  data  to  demonstrate  the  shift  in  m^^  detection  thresholds 

before  and  after  PHILTRE  was  applied  to  the  data.  Figures  24  and  25  show 
the  3— array  network  detection  rates  for  at  least  1 and  2 detections.  Table  VI 
shows  the  50  percent  detection  threshold  from  each  normal  curve.  Although 
insufficient  data  results  in  rather  wide  confidence  intervals  on  these  values, 
the  processing  results  generally  show  improvement  in  the  threshold  of  detec- 
tion. Considering  all  events,  there  is  a 50  percent  increase  in  the  number 
of  ALPA  detections,  a 61  percent  increase  in  the  number  of  LASA  detections, 
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Figure  21.  ALFA  — % detections  vs 
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Figure  24.  % detections  by  at  least  one  array  vs.  nr  for  a 3-array  network 


Figure  25.  % detections  by  at  least  two  arrays  vs.  iil  for  a 3-array  network 
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and  a 138  percent  increase  in  the  nunter  of  NORSAR  detections.  (Remember 

detections  were  attained  by  simply  bandpassing  NORSAR  data  ) 
There  are  94  percent  more  event,  detected  by  at  least  1 array,  44  percent  ' 
more  events  detected  by  at  least  2 arrays  and  67  percent  more  events  detected 
y arrays.  Since  the  A1P4  results  are  the  poorest  of  the  three  arrays, 
etectlon  rates  for  a 2-array  netuork  (LASA  and  NORSAR)  were  also  calculated. 
There  are  90  percent  more  events  detected  by  at  least  1 array  and  73  percent 
more  events  detected  by  2 arrays  in  this  case.  These  detection  results  show 
t at  the  thresholds  are  lowered,  perhaps  0.2-0. 3 magnitude  unit  (6  dB 
improvement).  Ihis  is  roughly  as  much  as  our  results  with  signals  burled  in 
noise  in  an  earlier  section  would  indicate.  We  note  that  spikes  in  the  AlPA 
beams  and  sudden  lumps  in  the  NORSAR  beams  produce  strange  effects  in  PHILTRE 
processing  and  were  responsible  for  many  spurious  signals  in  the  filtered 
output  which  hindered  detection  of  the  signal  of  Interest.  In  one  case 
PHILTRE  was  applied  to  an  ALPA  beam  of  a detected  LR  signal  that  contained 
several  spikes;  the  filtered  output  showed  no  detection.  The  quality  of  the 
input  data  can  thus  effect  the  ability  of  PHILTRE  to  retrieve  undetected  LR 
phases;  we  feel  that  the  results  could  have  been  significantly  better 
especially  for  ALPA.  if  the  data  input  to  PHILTRE  had  been  good  in  every  case. 
Latiert  et  al.  (1973)  report  somewhat  better  results  in  a similar  test  of  a 
PHILTRE-type  processor  on  LPE  data.  Lane  (1973)  reports  that  the  detection 
thresholds  for  beam  data  are  unchanged  after  PHILTRE  processing;  however,  that 

study  used  only  ALPA  beams,  with  which  our  results  were  the  poorest  of  the 
three  arrays. 


Lambert,  D. , Prahl,  S.,  and  Strauss,  A.,  1973,  Evaluation  of  the  noise 

detection  and  discrimination  capabilities  of  the 
--ions  and  L VLPE  NetwoS! 

three-component  filter,  Texas 
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CONCLUSIONS 


In  this  report  PHILTRE  showed  a mean  processor  gain  of  roughly  14  dB 
for  the  LR  phases  at  the  visual  detection  threshold  and  about  17  dB  at  5 dB 
above  the  visual  detection  threshold.  The  filter  should  be  able  t'  recover 
the  LR  phase  down  to  a S/N  input  of  about  -6  dB,  which  is  equivalent  to 
other  processors.  Further  research  may  show  that  the  PHILTRE  gain  is 
additive  to  that  from  match  filtering  or  maximum  likelihood  processing. 

Although  PHILTRE  deteriorates  rapidly  for  S/N  ratios  below  the  visual  detec- 
tion threshold,  its  significant  gain  at  the  visual  threshold  should  make  many 
marginal  detection  decisions  positive  near  the  threshold. 

PHILTRE  is  able  to  separate  two  signals  if  their  azimuthal  difference 
is  greater  than  about  60  degrees  and  if  at  the  same  time  the  amplitude  of 
the  second  signal  is  at  least  20  percent  o-'  the  amplitude  of  the  first  signal. 
Separation  is  possible  for  smaller  azimuthal  difference  only  if  the  signals 
have  reasonably  long  non-overlapping  parts. 

The  processor  PHILTRE  may  be  useful  in  recovering  from  ALPA,  LASA,  and 
NORSAR  beams  LR  phases  which  are  visually  undetectable  from  events  of  low 
body-wave  magnitude.  In  this  study  there  was  a 94  percent  Increase  in  the 
number  of  events  detected  by  at  least  one  of  the  three  arrays  and  a lowering 
of  the  50Z  threshold  by  0.2-0. 3 magnitude  unit. 

PHILTRE  should  be  applied  routinely  to  undetected  signals,  for  it 
should  recover  a significant  proportion  of  those  signals  which  are  below  the 
noise  level  or  which  are  interfered  with.  It  is  at  least  as  important  as 
bandpass  filtering  and  match  filtering  in  long-period  signal  detection  because 
PHILTRE  promises  at  least  as  much  gain  in  S/N  ratio  and  is  easy  to  apply,  although 
it  requires  more  computation. 
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